Translational Medicine @ UniSa, - ISSN 2239-9747 



2013, 7(4): 18-22 



Pharmacokinetics of Remifentanil: 
a three -compartmental modeling approach 

Sara Cascone 1 , Gaetano Lamberti* 1 , Giuseppe Titomanlio 1 , Ornella Piazza 2 

1 2 

Dipartimento di Ingegneria Industrial, Dipartimento di Medicina e Chirurgia, 
Universita di Salerno, Fisciano, Italy 



glamberti@unis a.it, www.minerva.unis a.it 



Abstract - Remifentanil is a new opioid derivative drug 
characterized by a fast onset and by a short time of action, 
since it is rapidly degraded by esterases in blood and other 
tissues. Its pharmacokinetic and pharmacodynamics 
properties make remifentanil a very interesting molecule in 
the field of anesthesia. However a complete and versatile 
pharmacokinetic description of remifentanil still lacks. In 
this work a three -compartmental model has been developed 
to describe the pharmacokinetics of remifentanil both in the 
case in which it is administered by intravenous constant-rate 
infusion and by bolus injection. The model curves have been 
compared with experimental data published in scientific 
papers and the model parameters have been optimized to 
describe both ways of administration. The ad hoc model is 
adaptable and potentially useful for predictive purposes. 

Keywords - Remifentanil, pharmacokinetics, three - 
compartmental model. 



I. INTRODUCTION 

The topic of the development of new opioid 
anesthetic agents is mainly to increase potency, and 
reduce the cardiovascular toxicity. For this purposes, 
recently a new kind of opioid derivative drug, 
remifentanil, has been synthetized. Remifentanil is an 
ultra-short acting opioid and it is subjected to metabolism 
by esterases in blood and other tissues. In vivo studies 
demonstrated an extensive metabolism of this drug by 
ester hydrolysis [1]. The primary metabolic pathway 
experienced by remifentanil is the formation of a 
carboxylic acid metabolite (named GI90291) obtained by 
de-esterification. The chemical structures of remifentanil 
and its primary metabolite are shown in Fig 1. It has been 
demonstrated on animals [1] that the pharmacodynamics 
of remifentanil is very similar to the other opioid drugs, 
this fact, combined with the reduced effects on the 
cardiovascular system makes the use of remifentanil very 
attractive in anesthesia. Remifentanil is generally 
administrated by the intravenous route. Because of the 
very short half-life of the drug, usually a bolus injection is 
administered to raise the blood concentration 
immediately, then a slower intravenous infusion is used to 
maintain the effective plasma concentration. Even if it is 
recommended to infuse remifentanil only during general 
anaesthesia procedures, the single/repeated bolus 
injections could be used in clinical situation in which a 
brief period of intense analgesia is required and the set - 
up of a continuous infusion pump is difficult (i.e. painful 
diagnostic and therapeutic procedures outside the 



operating theater). For this reason, it is particularly 
interesting to model what happen in the plasma 
concentration of remifentanil after the bolus injection or 
continuous infusion administration. Different kinds of 
models have been proposed, the simplest of which is the 
compartmental one [2], alternatively, the physiologically 
based approach [3] is more complex and potentially more 
exhaustive. 




Remifentanil GI90291 



Fig 1. Chemical structures of remifentanil and its metabolite. 

The aim of the present work is to develop and 
validate a new simple model using the conpartmental 
modeling approach to evaluate the remifentanil plasma 
concentration in case of bolus injection and continuous 
infusion. 

II. MODELING 

In the three-compart mental modeling, three 
compartments describe the fate of a drug once 
administered: the central compartment, which represents 
the plasma; the highly perfused compartment, which 
represents the organs and tissues highly perfused by the 
blood; and the scarcely perfused compartment, which 
represents the organs and tissues scarcely perfused by 
blood. A schematic of the model is shown in Fig 2. 
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Fig 2. A schematic of the three-compartmental model. 
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To model the pharmacokinetic of rernifentanil, the 
case of intravenous infusion has been studied. Thus, an 
amount of drug has been evaluated as inlet in the central 
compartment. The processes which cause the variation of 
the plasma concentration are: the absorption, the 
distribution, and the excretion of the drug. These 
phenomena have to be taken into account in the modeling. 
I(t) is defined as the function which describes the drug 
introduction by intravenous infusion (which could be an 
infusion at constant rate of administration or a bolus) in 
the central compartment. The drug concentration in the 
compartments could be evaluated solving the mass 
balance in the compartments, which could be written as: 



and bolus. Defining an error between the experimental 
data and the model evaluation: 



Vy 



Central compartment: 

" L = -Ch ■ C, + *21 -V 2 ■ C 2 + * 3 1 • C 3 ■ V 3 + 



dt 



V 2 



V3 
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Highly perfused compartment: 
dC 2 



dt 



= h 2 • C \ ■ V l - k 2 T C 2 ■ V 2 ~ Cl 2 ■ C 2 



Scarcely perfused compartment: 
dC 3 



dt 



■ = k l3 -Ci-Vi-k 3V c 3 -v 3 -a 3 -c 3 



In which, Cp, C 2 , and C3 are, respectively, the drug 
concentrations of the central, highly perfused, and 
scarcely perfused compartments. V7, V 2 , and V3 are, 
respectively, the volumes of the central, highly perfused, 
and scarcely perfused compartments. C/7, Cl 2 , and CI3 are, 
respectively, the clearances (rates of drug elimination) of 
the central, highly perfused, and scarcely perfused 
compartments. k 12 and k 21 are the transport coefficients 
between the central and the highly perfused 
compartments; k 13 and k 31 are the transport coefficients 
between the central and the scarcely perfused 
compartments. Finally, k 10 is the kinetic constant of drug 
elimination from the central compartment. The kinetics of 
elimination and transport between the compartments have 
been considered first order kinetics. These equations have 
to be solved coupled with their initial conditions: 



@t = 0,C P =0 

@t = 0,C 2 = 0 
@t = 0,C 3 = 0 



The three equations are inter-dependent, thus, they 
have to be solved simultaneously to evaluate the drug 
concentration in all the compartments. 

Once identified the transport phenomena which take 
place in the compartments and defined the differential 
equations able to solve the mass balances in the 
compartments, the value of the parameters has to be 
evaluated. To evaluate the parameters value, the model 
has been used to fit the experimental data taken in 
literature [4, 5], which refer both to intravenous infusion 



logQ 

In which n is the number of experimental data for 
each experiment, C H —C Bn (^ i ,p) is the difference 

between the experimental plas ma concentration and the 
model prediction at the time t i9 g_ being the vector of the 
model parameters. Minimizing this function it is possible 
to find the values of the parameters which better approach 
the experimental data. 

The set of three ODEs constituting the model have 
been solved by a code developed using Matlab, and the 
same software was used to find the value ofg_ minimizing 
the function e(£). 

m. RESULTS 

The model simulations obtained are shown in Fig 3 
and compared with the experimental data [5] in the case 
of intravenous constant-rate infusion. Each curve has been 
obtained as average values due to the administration to 
two subjects. During the 20 minutes infusion, a total of 14 
blood samples were taken. After stopping the infusion, 16 
blood samples were taken, up to 240 min after the stop. 
Therefore, each history was described by 30 sample data. 

As could be seen from the graphs, the model 
reproduces satisfactorily the experimental data. In 
particular, the plasma concentration following the 
administration of a dose of 1 ug-kg^-min" 1 is well 
approximated, however, the concentration of higher doses 
(4 and 8 ug-kg^-min" 1 ) are not well approximated for 
periods longer than 90 minutes, this is probably due to the 
fact that the supposed elimination kinetic could be still 
optimized. 

Furthermore, the same model has been used to 
reproduce the rernifentanil plasma concentrations in the 
case of bolus injection and compared with the 
experimental data [4] in Fig 4. In this case the 
administration has been evaluated as a fast infusion 
(bolus), thus the plasma concentration immediately rise to 
a high value. Each curve has been obtained as the average 
value over six patients (three men and three women). 
Over 360 minutes, 21 blood samples were collected and 
assayed for rernifentanil. Once again, the model curves 
satisfactorily approximate the experimental data. 
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Fig 3. Comparison between the experimental plasma concentration value 

[5] and the model curves in the case of intravenous constant-rate 
infusion with an infusion time of 20 minutes, ai) plasma concentration 
after a dose of 1 ug-kg" 1 -min" 1 ; sli) plasma concentration after a dose of 
4 ug-kg" 1 -min" 1 ; a3) plasma concentration after a dose of 8 ug-kg" 1 -min -1 . 



100 



0 10 20 30 40 50 60 70 80 90 100 

I.I.I L_ 



0.01 



b 2 > 




Dose = 2 MQ'kg 

■ experimental data 
model curve 



100 



10 



r1 



,0.1 



0.01 



10 20 30 40 50 60 70 80 90 100 
Time [min] 



10 20 30 40 50 60 70 80 90 100 
i i . i i i . i . i i i i i i i i i_ 



100- 



E 

Ui 

Ss 10- 
c 
o 



Dose = 5 jig- kg' 1 

■ experimental data 
model curve 



0.1 , 



0.01 



100 



10 



r1 



,0.1 



0.01 



0 10 20 30 40 50 60 70 80 90 100 
Time [min] 

Fig 4. Comparison between the experimental plasma concentration value 
[4] and the model curves in the case of fast intravenous infusion (bolus). 

bi) plasma concentration after a dose of 2 ug-kg" 1 ; 

hi) plasma concentration after a dose of 5 ug-kg 1 ; 

TABLE I. Values and dimensions of the three-compartmental model 
parameters. 



Parameter 


Optimized value 


Parameter 


Optimized value 


Vi 


7.88 mL 


fcio 


0.172 min" 1 


v 2 


23.9 mL 


k\i 


0.373 min" 1 


v 3 


13.8 mL 


ki\ 


0.103 min" 1 


Cli 


2.08 mL-min" 1 




0.0367 min" 1 


Ch 


0.828 mL-min" 1 




0.0124 min" 1 


Ch 


0.0784 mL-min" 1 







The values of the model parameters obtained after the 
optimization routine, are shown in Table 1. The model 
developed has been used to evaluate the plasma 
concentration both in the case of intravenous constant -rate 
infusion and intravenous bolus, which has been 
reproduced simulating a very fast infusion in the central 
compartment. This is a remarkable improvement to the 
compart mental modeling: in fact, once the model 
parameters have been evaluated for a certain 
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administration, the model is able to predict the drug 
plasma concentration varying not only the dose, but also 
the infusion rate of the drug. 
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Fig 4 (continue). Comparison between the experimental plasma 
concentration value [4] and the model curves in the case of fast 
intravenous infusion (bolus). 

b3) plasma concentration after a dose of 15 ug-kg 4 ; 

b0 plasma concentration after a dose of 30 ug-kg 1 . 

IV. CONCLUSIONS 

In this work a three -compart mental model has been 
developed to reproduce the evolution of remifentanil 
plasma concentrations after intravenous constant-rate 
infusion and intravenous bolus. The main phenomena of 
absorption, distribution, and metabolism have been 
identified and the mass balances for the three 
compartments have been written. The model has been 
then used to reproduce plasma concentrations taken from 
literature and the best values of the model parameters 
have been found minimizing the error between model 
curves and experimental data. 

Several studies have been conducted to develop a 
model which is able to reproduce the remifentanil 
pharmacokinetics. The aim of these studies is to compare 



the measured pharmacokinetic features of remifentanil 
after an intravenous infusion to the model prediction [6]. 
In particular, the compartmental analysis has been 
extensively used and compared with the experimental 
data, taken after intravenous infusion [5] or bolus [4]. 
Blood concentration and time data after a computer- 
controlled infusion of remifentanil could be analyzed by 
nonlinear regression using the NONMEM program 
(University of California) [7] which may produce 
predicted and individually predicted values (post hoc 
Bayesian estimates). The initial two-stage analysis 
comparing one-, two-, and three-compart mental models 
found that the two -compartmental model shows the best 
fit to the experimental data. This result was also 
confirmed by a population analysis. A more complex 
analysis of both pharmacokinetics and pharmacodynamics 
of remifentanil has been also approached [8]. The 
pharmacokinetic/pharmacodynamic relationship has been 
evaluated using non -linear regression analysis. The 
pharmacokinetics have been described using a one- 
compartment intravenous infusion model. Moreover, the 
pharmacodynamics have been fitted using inhibitory 
model. A statistical evaluation of the goodness of the 
models has been carried out following the Akaike 
Information Criterion [9]. According to this analysis, the 
decrease of the SSE using a model with a large number of 
parameters is useful if and only if it overcomes the 
increase in the number of parameters with respect to the 
use of a model with a limited number of parameters. In 
this case, the simple model has shown the best overall 
fitting results. 

Nevertheless, once the value of the parameters has 
been evaluated, our simple model was able to describe the 
remifentanil concentration on blood for different ways of 
administration. This is a remarkable improvement to the 
conpartimental modelling: in fact once the model 
parameters have been evaluated for a certain kind of 
administration, the model is able to predict the drug 
plasma concentration varying not only the dose but also 
the infusion rate of the drug. This feature makes the 
model more versatile than the other available in literature 
and very useful for predictive purposes. 
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